Bernoulli 14(3), 2008, 623-636 
DOI: 10.3150/08-BEJ128 



Asymptotic normality and consistency of a 
two-stage generalized least squares estimator 
in the growth curve model 

JIANHUA HU^ and GUOHUA YAN^ 

^Department of Statistics, Shanghai University of Finance and Economics, 777 Guo Ding Road, 
Shanghai 200433, P.R. China. E-mail: frank.jianhuahu@gmail.com 

^Department of Statistics, University of British Columbia, 333-6356 Agricultural Road, Vancou- 
ver, BC, Canada, V6T 1Z2. E-mail: guohuayan@gmail.com 

Let Y = X0Z' + f be the growth curve model with £ distributed with mean and covariance 
I„ (g) S, where 0, S are unknown matrices of parameters and X, Z are known matrices. For 
the estimable parametric transformation of the form 7 = C0D' with given C and D, the two- 
stage generalized least-squares estimator 7(Y) defined in (7) converges in probability to 7 as 
the sample size n tends to infinity and, further, y'n[7(Y) — 7] converges in distribution to the 
multivariate normal distribution A/'(0, (CR^^C) (g) (D(Z'S~^Z)"^D')) under the condition that 
lim„^oo X'X/n = R for some positive definite matrix R. Moreover, the unbiased and invariant 
quadratic estimator S(Y) defined in (6) is also proved to be consistent with the second-order 
parameter matrix S. 

Keywords: asymptotic normality; consistent estimator; estimation; generalized least-squares 
estimator; growth curve model 

1. Introduction 

The growth curve model is defined as 

Y = X0Z' + £, £~g(0,I„(8)S), (1) 

where Y is an 71 x p matrix of observations, X and Z are known n x to (n > m) and p x g 
[p > q) full-rank design matrices, respectively, is an unknown mx q matrix, called the 
first-order parameter matrix, and S is an unknown positive definite matrix of order p, 
called the second-order parameter matrix. £ follows a general continuous distribution Q 
with mean matrix and Kronecker product structure covariance matrix I„ (g) S. 

Model (1) was proposed by PotthofF and Roy [11] under the normality assumption of 
the error matrix £. Since then, parameter estimation, hypothesis testing and prediction 
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of future values have been investigated by numerous researchers, generating a substantial 
amount of literature concerning the model. 

In what follows, we give a brief review of the literature on large sample properties for 
the growth curve model, a particular kind of multivariate regression model. Chakravorti 
[2] presented the asymptotic properties of the maximum likelihood estimators. Zezula 
[15] investigated the asymptotic properties of the growth curve model with covariance 
components. Gong [4] gave the asymptotic distribution of the likelihood ratio statistic for 
testing sphericity. Bischoff [1] considered some asymptotic optimal tests for some growth 
curve models under non-normal error structure. However, no work has been done on the 
asymptotic normality and consistency of two-stage generalized least-squares estimators 
of the first-order parameter matrix for the growth curve model (1). 

In this paper, we shall investigate the consistency and asymptotic normality of a two- 
stage generalized least-squares estimator 7(Y) for the estimable parametric transfor- 
mation of the form 7 = C0D' with respect to the first-order parameter matrix 0. In 
addition, we shall demonstrate the consistency of a known quadratic covariance estimator 
S(Y) with the second-order parameter matrix S (see Zezula [14]). 

Readers are referred to Eicker [3], Theil [13] and Nussbaum [10] for results on the large 
sample properties of the least-squares estimators for ordinary univariate and multivariate 
regression models. 

This paper is divided into four sections. Some preliminaries are presented in Section 
2. In particular, for the estimable parametric transformation of the form 7 = C0D', a 
two-stage generalized least-squares estimator 7(Y) is defined in (7). The consistency of 
the estimator 7(Y) and the consistency of the known quadratic estimator S(Y) defined 
in (6) are investigated in Section 3. Finally, in Section 4, the asymptotic normality of the 
two-stage generalized least-squares estimator 7(Y) is obtained under a certain condition. 

2. Preliminaries 

Throughout this paper, the following notation is used. Let ^nxp denote the set of all 
n X p matrices. Let A' denote the transpose of the matrix A. Let tr(A) denote the 
trace of the matrix A. Let I„ denote the identity matrix of order n. For a sequence of 
numbers {a„} and a sequence of numbers we say that a„ = 0{bn) if there is a 

constant c such that limsup \a„/b„\ < c; we say that a„ = o{b„) if lima„/&„ = 0. For an 
n X p matrix Y, we write Y = [yi,y2, . . • ,yri]', Yi G 3?^, where 3?^ is the p-dimensional 
real space, and vec(Y') denotes the np-dimcnsional vector [y'i,y2, ■ . . ,y^]'- Here, the vec 
operator transforms a matrix into a vector by stacking the columns of the matrix one 
under another. Y ~ ^/(M, I„ (g) S) means that Y follows a general continuous distribution 
Q with E(Y) = M and that I„ S is the covariance matrix of the vector vcc(Y'); see 
Muirhcad [9], Section 3.1. The Kronecker product A® B of matrices A and B is defined 
to be A (g) B = (fly-B). We then have vec(ABC) = (C ® A) vec(B). Let A+ denote the 
Moore-Penrose inverse of A and Px = X(X'X)^X' be the projection onto the column 
space ^(X) of a matrix X along the orthogonal complement '^(X.)^ of '^^(K). 
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Given A S ^nxp and B G ^pxs, a linear parametric function B'/3 is called estimable 
with respect to A if there exists some T G ^nxs such that E(T'A/3) = B'/3 for all 
(3 gW; see Hu and Shi [6] for a more detailed description. 

Note that the first-order parameter in model (1) is defined before a design is planned 
and observation Y is obtained. Thus, the rows of the design matrix X in model (1) are 
added one after another and the term Z in model (1) docs not depend on the sample 
size n; see the example in Potthoff and Roy [11]. So, we shall only consider the case of 
full-rank matrices X and Z. 

As discussed in Potthoff and Roy [11], hypotheses of the form C0D' = under model 

(1) are usually considered, where C G ^sxm and D G ^txq- Thus, in this paper, we 
shall consider the estimator of the parametric transformation 7 = C0D' of with 
given matrices C S ^sxm and D £ ^txq- 

We shall begin by reviewing the case of a known second-order parameter matrix S, 
say Sq. According to the theory of least squares (sec, e.g., Rao [12], 4a. 2), the normal 
equations of model (1) arc X'X0Z'S(7^Z = X'YSo ^Z. The least-squares estimator ©o 
of is given by 

©0 = (X'X)-iX'YSo^Z(Z'So^Z)-\ (2) 

Since 

(X'X)-iX' = (X'X)-iX'Px 

and 

Z(Z'SoiZ)-iZ' = (PzSoiPz)+, (3) 

(2) can be written as 

©0 = (X'X)-iX'PxYSo-i(PzSo 'Pz)+Z(Z'Z)-i. (4) 

Let 

7o = C©oD'. (5) 

The mean and covariance of 7q are, respectively, C©D' and (C(X'X)~-'^C') ® 
(D(Z'SoiZ)-iD'). 

In addition, it follows from Rao [12], 4a. 2, that 7 = C©D', for any matrices C S ^sxm 
and D e ^txqi is an estimable parametric transformation if matrices X and Z are of 
full rank. So, 7q defined in (5) is said to be a least-squares estimator of the estimable 
parametric transformation 7 = C©D'. It is easily derived from 4a. 2 of Rao [12] that 7q 
is the best linear unbiased estimator (BLUE) of 7. 

Now, we shall focus our attention on the case of an unknown S. 

Let 

S(Y)=Y'WY, L^(i_p^.). (6) 

n — rank(Xj 

It is well known that I](Y)^^ is positive definite with probability 1 (see the proof of 
Muirhead [9], Theorem 3.1.4). Zezula [14], Theorem 3.7, tells us that S(Y) is a uniformly 
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minimum variance unbiased invariant estimator of S under the assumption of normality. 
This estimator S(Y) is often used to find the first-stage estimator; see, for example, 
Zezula [16]. We shall also take the estimator as the first-stage estimator in our following 
discussion. 

In (5), an unbiased least-squares estimator of 7 is given when S is known. However, 
when S is unknown, if we write = (X'X)-iX'YS~^Z(Z'S~^Z)"\ then the statistic 
7 = C0D' depends on S. In this case, we shall use a method called two-stage estimation 
to find an estimator, which is denoted by 7(Y): first, based on data Y, we find a first- 
stage estimator S of S; second, replace the unknown S with the first-stage estimator S 
and then find through the normal equations of model (1). 

We take S(Y) in (6) as the first-stage estimator S. Replacing S in (4) with S(Y), 
(5) can be expressed as 

7(Y) = C(X'X)-1x'Y£"\y)Z(Z'£"\y)Z)-1D'. (7) 

Let 

H(Y) = S"'(Y)(PzS"'(Y)Pz)+. (8) 

Then, by (3), (7) can be rewritten as 

7(Y) = C(X'X)-iX'YH(Y)Z(Z'Z)-iD'. (9) 

The estimator 7(Y) is said to be a two-stage generalized least-squares estimator of the 
estimable parametric transformation 7 = C0D'. 

In the special case of C and D being identity matrices, the estimable parametric 
transformation 7 is the first-order parameter matrix 0. By (9) or (4), we have 

0(Y) = (X'X)-iX'YH(Y)Z(Z'Z)~\ (10) 

The following lemma concerns the unbiasedness of the estimator 7(Y) under the as- 
sumption that S is symmetric about the origin. 

Lemma 2.1. Assume that the distribution of £ is symmetric about the origin. Then 
the statistic 7(Y) defined in (9) is an unbiased estimator of the estimable parametric 
transformation 7. 

Proof. Since S(Y) = = S(— £), 7(Y) can be expressed as 

7(Y) = C{X.'X.)-^X.'X&Z'f\e)Z{Z"E~\s)Z)-'^T)' 

+ c(x.'xy'^x.'st^\s)z{z't'\s)zy^T)'. 

Let 

M(£') = C{X'X.)-^X.'£t^\e)Z{Z't;^\s)Z)-^'D' . 
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Then M{-£) = -M{£) and hence E(M(5)) = 0. Thus, E(7(Y)) = C0D'. This com- 
pletes the proof. □ 



3. Consistency 

Since Y is associated with sample size n, we shall use Y„ to replace Y in (9) and then 
investigate the consistency of the estimator S(Y„), as well as the consistency of the 
related estimator 7(Y), as the sample size n tends to infinity. Note that X and £ are 
also associated with the sample size n. 

Recall that an estimator of S of the form Y^W*Y„ is unbiased and invariant if and 
only if tr(W*) = 1 and W*X = 0; see Zezula [14]. Hence, the statistic S(Y„) = Y^WY^, 
defined in (6) is an unbiased and invariant estimator of S without the assumption of 
normality. Moreover, under the assumption of normality, the estimator S(Y„) follows a 
Wishart distribution; see Hu [5]. 

Now, we shall investigate the consistency property of the estimator S(Y„). 

Theorem 3.1. For model (1), the statistic S(Y„) defined in (6) is a consistent estimator 
of the second-order parameter matrix S. 

Proof. Since Y;WY„ (Y„ - X0Z')'W(Y„ - X0Z'), in the following discussion we 
can assume without loss of generality that X0Z' = 0. So, by (6), 



S(Y„) = I - > ' £i£[ -S'J^xS I , (11) 



where 5 = (^i, £2, ■ • ■ , £n)' ^(0, I„ S). 

Note that {£i£i)f^i is a random sample from a population with mean E(£;£^J) = S. 
According to Kolmogorov's strong law of large numbers (see Rao [12], 2c. 3 (iv)), 

1 " 

— 7 5; 5; converges almost surely to S. (12) 
n ^-^ 

1=1 

Letting e > 0, by Chebyshev's inequahty and the fact that E(Y'WY) = tr(W)I] + 
E(Y)'W X E(Y), we have 



^ tr(I„ tr(S)Px) = ^ tr(Pjf ) tr(S). 









i 
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Since tr(Px) = rank(X) is a constant, P(||-^Px^^|l > e) tends to as the sample size 
n tends to infinity. So. 

—=Px£ converges in probability to 0. (13) 

Since convergence almost surely implies convergence in probability, by (12) and (13), we 
obtain from (11) that S(Y„) converges to S in probability. This completes the proof. □ 

Now, wc focus our attention on the consistency of the estimator 7(Y,i). We first prove 
the following lemma. 



Lemma 3.2. H(Y„) converges in probability to H, where H(Y„) is defined in (8) and 
H = S-i(P2S"iPz)+. 



Proof. Note that the function A to A+ is not continuous. Since S (Y^) is positive 
definite with probability 1, by Lehmann [7], Lemma 5.3.2, and Theorem 3.1, we have 



S (Y„) converges in probability to S ^ . 



(14) 



Write 



Pz = OAO', Q„ = 0'S(Y„)0, 

where O is a p x p orthogonal matrix, A = diag[0, Ig] with q = rank(Z) and 



Q-^^O't (Y„)0 = 



Gii(Y„) Gi2(Y„) 
G2i(Y„) G22(Y„) 



— [G.y(Y„)]2x2 



with G22(Y„) a qx q random matrix. By (14), for any i,j — 1,2, Gij(Y„) converges in 
probability to Gij. Note that (PzCPz)+ = Pz(PzCPz)+Pz- Then 

H(Y„) = S"\y„)(PzS"\y„)Pz)+ = OQ,-iAO'(OAQ-iAO')+OAO' 
= 0(G„-(Y„))2x2AO'(OA(G,,(Y„))2x2AO')+OAQ' 



O 



o 



Gi2(Y„) 

G22(Y„) 

Gi2(Y„) 

G22(Y„) 



0'(0 diag[0, G22(Y„)]0')+OAO' 



diag[0,G2-2i(Y„)]AO' = O 



G*(Y„ 
I„ 



O', 



where G*(Y„) = Gi2(Y„)G^2'(Y„). 
Similarly, H can be decomposed as 

H = 



Gi2G22^ 



O'. 
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G*(Y„ 

la 



G12G22 
I„ 



converges in probability to 
namely, H(Y„) converges in probability to H. This completes the proof. 



□ 



Based on Lemma 3.2, we obtain the following consistency result for the estimator 
7(Y„). 



Theorem 3.3. At 



that 



lim -(X'X) =R, 

n— >cx) n 



(15) 



where R is a positive definite matrix. Then the statistic 'y{Yn) is a consistent estimator 
of the estimable parametric transformation 7 = C0D'. 

Proof. To prove that 7(Y„) is a consistent estimator of 7, by Slutsky's theorem (see 
Lehmann and Romano [8], Theorem 11.2.11), it suffices to show that 0(Y„) is a consis- 
tent estimator of 0. 

Replacing Y with X0Z' + £ in (10), wc decompose 0(Y„) as E„ + F„, where 



(X'X)"iX'X0Z'H(Y„)Z(Z'Z)-i 



and 



F„ = (X'X)-^X'£H(Y„)Z(Z'Z)-\ 



Note that A'A(A'A) B = B if B is estimable with respect to A. Since Z'P^ = Z', we 
obtain 

E„ = 0Z'S"'(Y„)(P2S"\y„)Pz)+Z(Z'Z)-i 

= 0Z'PzS"\y„)Pz(PzS"\y„)Pz)+Z(Z'Z)-i 
= 0Z'Z(Z'Z)-i =0 



and 



F„=n(X'X)- 



=X' 



Px5 H(Y„)Z(Z'Z)- 



By (15), Xy-^n are bounded. In fact, the elements of X'/y^ are at most of order n~^/^ 
(see the proof of Lemma 4.1 below). Then, by (13), (15) and Lemma 3.2, F„ converges in 
probability to 0. So, 0(Y„) converges in probability to 0. This completes the proof. □ 
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4. Asymptotic normality 

We investigated the consistency of the estimator 7(Y„) in Section 3. In this section, we 
shall investigate the asymptotic normality of y^[7(Y„) — 7]. 
First, we shall prove the following lemma. 

Lemma 4.1. Suppose that condition (15) holds. Let S = (X'X)^^X' = (si, S2, . . . , Sn)rnxn, 
where S; is the Ith column of (X'X')~^X'. Then, for any I G {1, 2, . . . , n}, the m elements 
of y/nsi are 0[n^^/'^) . 

Proof. Write V= ^^X' = [vi, V2, . . . , v„]. The transpose of v; is an m-clement row 
vector, 

11 1 \ 

zXii, —=Xi2, . . . , -^Xi,n , 



,-1 v'-" 



where X — [xyjnxm- By (15), W = ri^^X'X converges to a positive definite matrix R. 
So, the elements of VV — + V2V2 + • • • + v„v'j are bounded. We claim that for any 
I £ {1, 2, . . . , n}, the m elements of v; are all 0(rr^l^). 

If this is not true, we can assume, without loss of generality, that one element of v„ 
is OinP'^l"^) with p > 0. Then one element of v„vjj would be 0(ri^^~^). Hence, the 
corresponding clement in matrix VV = Viv'j^ + V2V2 + • • • + v„v'^ would be O(n^J'), 
which is not bounded. This contradicts condition (15). 

Note that 



(V^si, V^S2, . . . , V^s„) = x/^(X'X)-^X' = n(X'X)-i^X' 



n 



= 7l(X'X) ^[V1,V2, . . . , v„], 

namely, for / = 1, 2, . . . , n, ^/ns; = ?i(X'X)^^v/. Thus, for ^ = 1, 2, . . . , n, the m elements 
of -y/ns; are also 0{n^^l^\ This completes the proof. □ 

Now, we shall show the following important result on the asymptotic normality of 
V^[7(Y„)-7]. 

Theorem 4.2. Under the assumption of condition (15), y/n['y{'Yn) ^7] converges in 
distribution to Afsxt{0, {CR-^C) (g> (D(Z'S"^Z)-iD')). 

Proof. First, by (8), we rewrite 7 and 7(Y„) as 

7 = CSX0Z'H(Y„)KD' 

and 

7(Y„) = CSY„H(Y„)KD', 
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where K = Z(Z'Z)-i. So, 

7(Y„) - 7 = CSY„H(Y„)KD' - CSX0Z'H(Y„)KD' 

= CS(Y„-X0Z')H(Y„)KD' (16) 
= CS£H(Y„)KD' = CL„H(Y„)KD', 

where L.^^ = SS. Further, is expressed as 

n 

L„=^s,£;, (17) 

1=1 

where s; is the Ith column vector of S and £[ is the Ith row vector of the matrix £ with 
£~e(0,I„®S). 

Next, we shaU find the hmiting distribution of ^/n['y(Yn) — 7] through showing that 

y/nhn converges in distribution to A/'„ixp(0, (X) S). (18) 

Since {£'i}"^i are independent and identically distributed, for t g ^pxm, the charac- 
teristic function ^'n (t) of y^L'^ is given by 

n 

4'„(t) = E(cxp{itr(t'V^L;)}) =E(exp{itr(tV^L„)}) =[|$(\/^tsO, 

where $(•) is the characteristic function of 
Recall that for u in the neighborhood of 0, 

\n{l - u) = -u + f{u) with f{u)^^u^ + o{u^). (19) 

If we write p{u) = f(u)/u, then from (19), 

p{u)=o{u) asit— >0. (20) 

Also, 

$(x) = l-ix'Sx + .g(x) forxeSR™ and ,g(x) o(||xj|2) asx^O. (21) 
For e > 0, there exists S{e) > such that 

|g(x)| <e||x|p as 0< l|x|| <(5(£). (22) 
Therefore, by (19) and (21), the characteristic function of ^/ni^'n can be decomposed as 



4'„(t) = exp <^^ln($(V^tsO) 



. 1=1 
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(23) 



expi 



. 1=1 



--nsJt'Stsj + g{^/ntSl) + / -nsJt'Sts; - g{y/nts 



where 



^1=1 



a„ = ^ nsJt'Sts; = tr ( ^ nsJt'Ets; ] , 

1=1 

n 



and 



For an, we have 



By (15), 



1=1 



Vn = Y /(insJt'Sts/ - g{^/ntSl)). 

1=1 



an = tr(StnSS't') = tr(Stn(X'X)-H'). 



hm an = tr(R"H'St) (vcc(t))'(R~^ (g) S) vcc(t) 



(24) 



(25) 



For Pn, by Lemma 4.1 and the contmuity of ts;, for the S{e) > in (22), there is an 
integer N{e) > such that for n > N{e), 



< II V^ts/|| < (5(e) for alU = 1, 2, . . . , n. 
If we take n > N{e), then by (22) and (26), 

|5(V^ts;)| < llVntsjIl^e. 

So, 

71. 

|/3„ I < ^ II V^ts; II \ = e tr(ntSS't') = £ tr(t7i(X'X)-H'). 



(26) 



(27) 



1=1 



So, by (15), hmsup„^o^ |/3„| < £tr(tR t'). Since e > is arbitrary, we obtain 



hm (3n = 0. 



(28) 
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So, by (27), 



A; = i(Vn^ts/)'S(Vnts/) - .g(^A^ts^). 



\Xi\ < i(V^ts,)'S(V^ts,) + WV^tsife. 



(29) 



Take n > N{e). By Lemma 4.1, the continuity of ts; and (20), increasing N{e) if necessary, 
we may suppose that for all I, \p{Xi)\ < e. Since /(A;) =p{Xi)Xi, 



1=1 



So, by (29), 



or 



\Vn\<J2 

1-- 

n 

l'?n|<E 



■tr(V^sJt'Sti/ns;) + ||Vnts/| 



tr(t'Sts/s;) +tr(V?its;(Vnts/) £^) 



namely. 



\Vn\ < -tr(t'StriSS') +tr(tnSS't')e2. 



Note that nSS' = n(X'X)"i. Since e is arbitrary, by (15) and (30), 

lim rjn = 0. 

n — ^00 

By (25), (28) and (31), we obtain from (23) that 

lim «'„(t) =exp{-i(vec(t))'(R"i(8)S)vec(t)}. 



(30) 



(31) 



(32) 



So, by Levy's continuity theorem, y/nLn in (17) converges in distribution to the normal 
distribution A/'mxp(0,R~^ ® S), as was claimed in (18). 

Finally, by Lemma 3.2, (16), (18) and Muirhead [9], Theorem 1.2.6, we obtain that 

V^[7(Y„) - 7] converges in distribution to 7V;,xt(0, CR^^C ® (DK'H'SHKD')). 

Replacing H and K with S"^(PzS"^Pz)+ and Z(Z'Z)-i, respectively, we conclude 
that 

Vn[7(Y„) - 7] converges in distribution to 7Vsxt(0, (CR^^C) ® (DTD')), 

where T = (Z'Z)-iZ'(Pz5]-ip2)+Z(Z'Z)-i = (Z'^-^Z)"! (see (3)). Thus, the proof 
is complete. □ 
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Under model (1), hypotheses of the form 

H:7 = C0D' = O 

are usually considered; see PotthofF and Roy [11]. 

From Theorem 4.2 and Slutsky's theorem, the following corollary provides the asymp- 
totic behavior of ^/n7(Y„) under H. 

Corollary 4.3. Under the assumption of condition (15), if matrices C(X'X)^^C' and 
D(Z'S"\y„) Z)-1D' are non-singular, then the statistic (Cn(X'X)-iC')^^/^\/^7(Y„)(D(Z' x 
tr (Y„)Z)-iD')"^/^ under H converges in distribution to MsxtiO,T) . 

Remark 4-4 • Lemma 2.1 tells us that 7(Y) is an unbiased estimator of 7 under the 
assumption of £ being symmetric about the origin. In general, it is very difficult to 
obtain the covariance matrix of 7(Y), even under the assumption of normality. How- 
ever, under condition (15), Theorem 4.2 gives us an approximate covariance matrix 
(C(X'X)-iC') (g> (D(Z'S"\y)Z)-iD') of 7(Y) for large sample size n, without the 
assumption of normality. 

We now conclude this paper by discussing the example in PotthofF and Roy [11]. No 
assumption of normality is made in our discussion. 



Example 4-5. There are to groups of animals, with r animals in the jth group and each 
group being subjected to a different treatment. Animals in all groups are measured at the 
same p time points, ti,t2, ■ ■ ■ ,tp. The observations of different animals are independent, 
but the p observations on each animal are assumed to have a covariance matrix S. 

Based on the problem and our discussion, m remains constant, while r tends to infinity. 

For i = 1, 2, . . . , TO, the growth curve associated with the ith group is 



'jO + OiiX + 9i2X^ H 1- Oiq-iX' 



9-1 



Put 



where = [5ue'^,52ie'^, 
n = rm, 



X= (X1,X2,...,X„), 

,Sue'^]',er^il,l,...,iyeW 



e, 



7iO, (7^1,(7^2, 



5ij are the Kronecker symbols. 



@ = {e[,e'2,...M' 



and 



Z = 



t2 



tl 



+9-1 



. 1 tr, ti 
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The observation data matrix Y„ can be written as 

Y„ = X0Z' + £, 

where £ = {£1,62, . ■ . ,£nY with E\,Si^ ■ ■ ■ being independent and identically dis- 
tributed with mean and covariancc S. Then, by (14), 

R-i = lim ri(X'X)-i =toI. 

r — *-oo 

By (10), 

771 „ 1 „ 1 

0(Y„)--X'Y„£ (Y„)(Pz(S (Y„))Pz)+K. 

n 

For the estimable parametric transformation of the form 7 = C0D' with given C € 
^sxm and D g ^txg, the two-stage generalized least-squares estimator is given by 

7(Y„)-C0(Y„)D'. 

It follows from Theorem 4.2 that ^/?i[7(Y„) — 7] converges in distribution to the normal 
distribution A/;xt(0, (mCC) O (D(Z'S"iZ)-iD')). 

Moreover, if wc try to test that all m growth curves are equal, except possibly for the 
additive constant ^iOj then we take C to be a matrix whose last column contains all — I's 
and whose first (m— 1) columns constitute the identity matrix, and D to be a (g — 1) x q 
matrix whose first column contains all O's and whose last (g — 1) columns constitute the 
identity matrix, namely, taking 

C -lm-l](m_i)xm > D=[0 Ig-l](q_i)xg ' 

where Im-i ~ (I7 li ■ • ■ , 1)': and hypothesis Hp : C0D' = 0. Obviously, matrices 

C(X'X)-iC' and D(Z'S~\y„)Z)-iD' are non-singular. It follows from Corollary 4.3 

that statistic (Cn(X'X)-iC')"^/^V^7(Y„)(D(Z'S~\Y„)Z)-iD')"^/2 under Ho con- 
verges in distribution to A/'(m-i)x((j-i)(0,I). 
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